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Abstract 



Some dynamical properties present in a problem concerning the acceleration of particles in a wave packet are studied. The dynamics 
^"™^ of the model is described in terms of a two-dimensional area preserving map. We show that the phase space is mixed in the sense 
that there are regular and chaotic regions coexisting. We use a connection with the standard map in order to find the position of 
the first invariant spanning curve which borders the chaotic sea. We find that the position of the first invariant spanning curve 
increases as a power of the control parameter with the exponent 2/3. The standard deviation of the kinetic energy of an ensemble 



CO 



of initial conditions obeys a power law as a function of time, and saturates after some crossover. Scaling formalism is used in order 



to characterize the chaotic region close to the transition from integrability to nonintegrability and a relationship between the power 
- s I law exponents is derived. The formalism can be applied in many different systems with mixed phase space. Then, dissipation is 
introduced into the model and therefore the property of area preservation is broken, and consequently attractors are observed. We 
i— i show that after a small change of the dissipation, the chaotic attractor as well as its basin of attraction are destroyed, thus leading 
(i_J| the system to experience a boundary crisis. The transient after the crisis follows a power law with exponent -2. 
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During the last decades many theoretical studies of area- 
preserving maps have intensively been done [[j] IH. One of 
the most studied models is the standard map proposed by B. V. 
Chirikov in 1969. The model describes the motion of the 
kicked rotator. The standard map can also be applied in differ- 
ent fields of science including accelerator physics |5|], plasma 
physics 1 6] and solid state physics [7]. It has also been stud- 
ied in relation to problems of quantum mechanics and quantum 
chaos fslfi OH, statistical mechanics fill] and many others. 

The system we are considering in this paper is the dynam- 
ics of a particle under the action of electrostatic waves. Usually, 
the model is described using the formalism of discrete nonlin- 
ear maps. It is well known that the structure of the phase space 
in such systems depends on the combinations of both initial 
conditions and control parameters. Basically they are classi- 
fied in three different classes: (i) integrable, (ii) ergodic and 
(iii) mixed. In case (i), the phase space consists of invariant 
tori filling the entire phase space. In case (ii), the ergodic sys- 
tems have only one chaotic component, and the time evolution 
of a sing le initial condition is enough to fill the phase space 
1511 . Finally, the case (iii) is most common and an im- 
portant property in the mixed phase space is that chaotic seas 
are generally surrounding Kolmogoro v- Arnold-Mo ser (KAM) 
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islands which are confined by a set of invariant curves. These 
mixed type systems are subject of intense research in classical 
and quantum chaos d-[32]. The structure of the phase space 
of the model we are considering is mixed and we derive an- 
alytically the location of the invariant curves. We consider a 
connection with the standard map close to the transition from 
local to global chaos and obtain an effective control parame- 
ter as well as the location of the first invariant torus bordering 
the chaotic sea in the phase space. In this paper we study the 
scaling behaviour of the model close to the transition from in- 
tegrability to nonintegrability [l33ll in a similar way as in [[34ll . 
as well as the dissipative version, and, moreover, we analyze 
the boundary crisis. Our main concern is the acceleration of the 
particle J35l [3(1 [3T|] and the saturation of the velocity [l38ll . An 
extensive work using scaling arguments was published recently 
in D4fl . However, when the dissipative dynamics is taken into 
account the structure of the phase space is changed and attrac- 
tors, namely, attracting fixed points and chaotic attractor, might 
be observed. Then, increasing the strength of the dissipation 
the edge of the basins of attraction of the attracting fixed points 
and the edge of the basin of attraction of the chaotic attractor 
(which corresponds to the stable and unstable manifolds of a 
saddle fixed point) touch each other and as a result the chaotic 
attractor as well as its basin of attraction is immediately de- 
stroyed. Such an event is called boundary crisis ||39l|40|]. After 
the crisis, the particle experiences a chaotic transient in the cor- 
responding region where a chaotic attractor existed prior to the 
crisis. The transient is described by a power law with respect 
to the relative distance in the control parameter where the crisis 
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happens. 

The paper is organized as follows. In section [2] we con- 
struct the two-dimensional map that describes the dynamics of 
the system and discuss our numerical results. Conclusions are 
drawn in section [3] 

2. The model and numerical results 

We shall study the model of a point mass (m) charged (e) 
particle moving in an electric field wave packet, as introduced 
by G. Zaslavsky et al. 114 In . which in general can be written as 



e 
m 



^ E n cos(k n x - TU n t), 



(1) 



where E n is the amplitude of the n-th Fourier component of 
the electric field wave. They consider a wave packet with the 
broad spectrum, and furthermore assume E n = Eq for all n. 
Moreover, rn n « tuq + nAm and k n « & + nAk. The group 
velocity of the wave packet is thus v g = dm/dk « Am/Ak. If 
the particle's velocity v = x is much larger than v g , one can use 
the approximation m n ^ mo. Using all these assumptions we 
get 



x = — Eq cos^o* - tcJoO ^ cos(nAkx), 



(2) 



and therefore using the Fourier decomposition of the periodic 
Dirac delta function we have 



x = —LEq cos(kox - mot) ^ S(x - nL), 



(3) 



where we use Ak = 2n/L. Between the delta kicks we have 
free motion of the particle and therefore the system of the two 
first order ordinary differential equations for x and v emerging 
from © can be formulated exactly as a discrete map. Defining 
6 = kox - TUot and 77 = mv\v\/2 (kinetic energy), and denoting 
by index n their values just before the n-th kick, we can write 
this map in the following form 



P : 



r] n+ i = rj n + eE Lcos(6 n ) 

6 n+ i =6 n + k L sign(7/ w+ i) - m L ^m/2\r] n 



,(4) 



where sign(7/ n+1 ) = 1 if rj n+ i > or sign(7/ n+1 ) = -1 if 77^+1 < 
0. However, defining = Irilmm^l} and j3 = 0/2n - 1/4 as 
auxiliary variables, and introducing a dissipation parameter <r, 
the map P can be rewritten as 



Vn+l 



(1 -o-)cf) n + 7 sin(27r^) 



Pn+l - [fin ~ 



1 



27T\(P n+l \ 



- + ^sign(0 n+ i)] mod(l) 



,(5) 



where y - leEo/mur^L, cr e [0, 1] and S = koL/2n. It is im- 
portant to stress that 6 is just a shift of the phase and from now 
on we fix it as S = without loss of generality. Here, y is the 
control parameter which controls the transition from integrabil- 
ity (y = 0) to nonintegrability (7 ^ 0). As shown by Zaslavsky 
et al. ll4lll for large 7 » 1 we have adiabatic picture with in- 
termittent but largely regular behaviour. The determinant of the 
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Figure 1: (Color online) (a) Phase space generated from Eq. [5] the control 
parameter were 7 = 10 -5 ; (b) Convergence of the positive Lyapunov exponent. 



Jacobian matrix is det J = 1 - cr. So for the case of cr = the 
mapping is area preserving. Here we shall investigate some dy- 
namical properties close to the transition from integrability to 
nonintegrability assuming 7 « 1. 

Let us first consider the conservative case where cr = 0. Fig- 
ure \T\ (a) shows the phase space generated from Eq. © where 
the control parameter used is 7 = 10" 5 . As one can see the 
phase space is mixed in the sense there are regular regions, in- 
variant spanning curves (invariant tori) and KAM islands coex- 
isting with a chaotic sea. The Lyapunov exponent is an impor- 
tant tool that can be used to classify orbits as chaotic or not. As 
discussed in |42|], the Lyapunov exponents are defined as 



1 



lj = lim - In |A 7 -| , j 

n—>c*3 n 



1,2 



(6) 



where Aj are the eigenvalues of M = n?=i Jiifiu <Pd an d h is 
the Jacobian matrix evaluated over the orbit {fii, (pi). However, a 
direct implementation of a computational algorithm to evaluate 
Eq. © has a severe limitation to obtain M. Even in the limit of 
short n, the components of M can assume very different orders 
of magnitude for chaotic orbits and periodic attractors, yielding 
impracticable the implementation of the algorithm. In order to 
avoid such problem, we note that / can be written as / = ®T 



2 



where © is an orthogonal matrix and T is a right up triangu- 
lar matrix. Thus we rewrite M as M = J n J n -\ • • • Ji®\®\ l J\, 
where Ti = &[ l J\. A product of /2©i defines a new m a 
next step, it is easy to show that M = J n J n -\ . . . /3©2©2 1 ^2^i- 
The same procedure can be used to obtain T2 = ©2^2 a so 
on. Using this procedure, the problem is reduced to evaluate 
the diagonal elements of T t : T[ v T l ir Finally, the Lyapunov 
exponents are given by 



1 n 
lim - V In \T 

i=\ 



1,2 . 



(7) 



If at least one of the Aj is positive then the system is classified 
as chaotic. Figure [T] (b) shows the behaviour of the positive 
Lyapunov exponent. We have used 10 different initial condi- 
tions iterated in the large chaotic region up to 10 9 times. The 
average of the positive Lyapunov exponent for our ensemble is 
A = 1.131(8), where 0.008 corresponds to the standard devia- 
tion of the 10 samples. 

The structure of the phase space of the discrete dynamical 
system we are dealing with and defined by the map P [Eq. (0)] 
is illustrated in Fig. [U(a). We can make a connection between 
the standard map and the map P in order to find the position of 
the first invariant spanning curve bordering the chaotic region 
and also to characterize the transition from integrability to non- 
integrability as it has been done in II2[ |43ll . The standard map 
was first proposed by J. B. Taylor [|44ll in order to study the 
existence of the invariants of motion in magnetic traps. Later, 
B. V. Chirikov ||3|,|4|] proposed another way to obtain the map 
and since then it became clear that situations described by the 
Chirikov map occur in many physical systems. The Chirikov 
map is 



C : 



I n+i + sin (OJ 

o w+1 = a> n + I n+1 



(8) 



where K is the control parameter (kick parameter). K controls 
the transition from integrability, (K - 0), to nonintegrability, 
(K ± 0), and also controls the transition from local chaos K < 
K c , where there still exists a set of invariant spanning curves 
separating different regions in the phase space, to global chaos, 
K > K c , where all the global invariant curves are destroyed. 

The critical K is K c = 0.971635 The connection of this 

result with the present model © is as follows. Suppose that 
close to the first invariant spanning curve, cp n can be written as 



1>n = <P* + A0n 



(9) 



where 0* is a typical value along the invariant spanning curve 
and A(p n is a small perturbation of 0. After defining X n = 2nf3 n , 
the second equation in the map © can be written as 



X n+ ] — X n — 



1 



Vn+l 



1/2 



Using Eq. ©, we can rewrite Eq. (TTTTb as 

1-1/2 



X n +\ — X n ■ 



1 



|l/2 



1 + 



A0, 
|0* 



n+l 
1/2 



(10) 



(11) 
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Figure 2: (Color online) (a) Behaviour of the standard deviation as function 
of n for different values of the control parameter 7. (b) Their initial collapse 
after the transformation ny 1 . (c) Their collapse onto a single and universal plot. 
The different curves are coded by the symbols in the legend in the same order 
(top-down). 



Expanding Eq. (fTTb in Taylor series and taking into account 
only terms of first order, we can rewrite Eq. (ITU as 



X n+ \ — X n + 



11/2 



1 - 



2101 



(12) 



Using Eq. ©, the first equation of the map © can also be 
written as 



+ A(f) n+ i = 0* + A(p n + ysm(2np n ) . 



(13) 



Multiplying both sides of Eq. (TT31) by ^ and adding -pjiyi, 

(14) 



we define 
/» = 



1 + A(f>n 



|0*|i/2 2 |0*p/ 2 ' 
and re- write the map © as 

g . 1 In+i = In + 2|0^|3/2 s i n (^«) 
X n +i = X n + I n+ \ 



(15) 
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z=-0.657(3) 
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Figure 3: (Color online) (a) Plot of aj sat as function of the control parameter y. 
(b) Behaviour of the crossover number n x against y. 



Comparing the standard map [see Eq. (|8])] and the map S , we 
see that there is an effective control parameter K e ff which is 
given by 



K eff ~ 



2|0*| 3 / 2 ' 



(16) 



Since the transition from local to global chaos occurs at K e ff > 
0.971635 . . ., the location of the first invariant spanning curve 
is given by 



r 



\2/3 



• • (17) 
^2x0.971635.../ v J 

According to Eq. ([T71) , the position of the first invariant span- 
ning curve changes with exponent 2/3 as the control parameter 
y varies. We conclude that for a given y a large chaotic sea is 
bordered by | |<| 0* |. For example, for y = 10" 5 , as in Fig. 
[Ha), we see indeed that 0* = 2.9 x 10" 4 in perfect agreement 
with Eq. ([T7l) . Similarly as in the standard map in such regime 
I 2 grows diffusively, i.e. I 2 oc n 114311 , 

We now assume small values for the parameter y and study 
the average standard deviation of (kinetic energy), which is 
defined as 



coin, y) 



cp 2 (n,y)-cpi in,y), 



(18) 



cp(n, y) ■■ 



1 n 

i=0 



(19) 



We have iterated Eq. (IT8T> up to n = 10 8 for an ensemble of C = 
1000 different initial conditions. The variable 0o was fixed as 
0o = 10" 2 y andyfto was uniformly distributed on [0, 1]. Figure [2] 
shows the behaviour of coin, y) as function of n for five different 
values of the control parameter y, as labelled in the figure. It is 
easy to see in Fig. [2ja) two different kinds of behaviour. For 
short n, coin, y) grows according to a power law and suddenly 
it bends towards a regime of saturation for long enough values 
of n. The crossover from growth to the saturation is marked by 
a crossover iteration number n x , which is very well defined by 
the intersection of the acceleration and saturation straight lines. 
It must be emphasized that different values of the parameter y 
generate different behaviors for short n. However, applying the 
transformation n — > ny 2 all the curves start growing together for 
short n, as shown in Fig. EJb). Based on the behaviour shown 
in Fig. [2l we can suppose the following scaling hypotheses: 



1. When n <?c n x , coin, y) grows according to 
coiny 2 ,y) oc iny 2 "f , 



(20) 



where the exponent (5 is the acceleration exponent; 

2. For a long number of iteration, n » n x , coin, y) approaches 
a regime of saturation which is described by 

oJ S atiny 2 ,y)ocy\ (21) 

where the exponent a is the saturation exponent; 

3. The crossover iteration number that marks the change 
from growth to the saturation is written as 



n x ocy 1 , 

where z is the crossover exponent. 



(22) 



After considering these three initial suppositions, we describe 
coin, y) in terms of a generalized homogeneous function of the 
type 



coiny 2 , y) = Tcoir a ny 2 , r b y) , 



(23) 



where r is the scaling factor, a and b are scaling exponents that 
in principle must be related to a, J3 and z. Since r is a scaling 
factor we can choose r a ny 2 = 1, or r = iny 2 )~ lla . Thus, Eq. 
(l23t is rewritten as 



coiny 2 , y) = iny 2 y l,a co x \iny 2 y b,a y\ , 



(24) 



where the function co\[iny 2 )~ b/a y] = co[l,iny 2 )~ b/a y] is as- 
sumed to be constant for n <?c n x . Comparing Eqs. (l20b and 
(|24L gives us ft = -I /a. Choosing now r b y = 1, we have 
T = y~ l l b and Eq. (f23l) is given by 



coiny 2 , y) = y- l/b co 2 iy~ a/b ny 2 ) , 



(25) 



where co^iy alb ny 2 ) = coiy alb ny 2 , 1). It is also assumed as 
constant for n » n x . Comparison of Eqs. (I2TT) and (l25t gives us 




Figure 4: (Color online) Characterization of a boundary crisis. The control parameters used were y = 10 2 and: (a) <x = 0.15 (immediately before the crisis); (b) 
cr = 0.143 (immediately after the crisis). 



a = -l/b. Given the two different expressions of the saturation 
value (jixj 1 ^ oc y a , we obtain 

-2 



n x oc yP 



(26) 



Thus comparing Eq. (l26l) and Eq. (l22k the crossover exponent 
is z = a IP - 2. Note that the scaling exponents are determined 
if the critical exponents ft and a were numerically obtained. The 
exponent J3 is obtained from a power law fitting for co(ny 2 ,y) 
curves for the parameter y e [10~ 9 , 10" 5 ] for short iteration 
number, n « n x . Thus, the average of these values gives us/3 = 

0. 499(4) = 1/2. Figure [3] shows the behaviour of (a) a) sat vs.y 
and (b) n x vs.y. Applying power law fitting in the figure, we 
obtain a = 0.6736(8) = 2/3 and z = -0.657(3). We can also 
obtain the exponent z evaluating z = a//3 -2 and the previous 
values of both a and J3. We found that z = -0.651(9). This 
result indeed agrees with our numerical data. This set of critical 
exponents makes this transition to belongs to the same class 
of universality as the phase transition observed in the problem 
of a classical particle confined inside an infinitely deep box of 
potential that has an oscillating square well 114611 or time varying 
barrier in the middle fl47l1 . 

Finally, in order to confirm our initial hypotheses and, since 
the values of the critical exponents a, J3 and z are now already 
known, we can use the scaling laws. Fig. [2jc) shows five differ- 
ent curves of oj generated from different values of the control 
parameter y overlapped onto a single universal plot. 

Let us now present our results for the dissipative dynamics 

1. e. cr ^ 0. It is well known that when dissipation is considered, 
the mixed structure of the phase space is changed and ellipti- 
cal fixed points can in principle turn into attracting fixed points 
(sinks) and chaotic seas might be replaced by chaotic attractors. 
Figure|4]shows the corresponding stable and unstable manifolds 



for two different saddle fixed points which are obtained by solv- 
ing (p n+ i = (p n and p n+ \ - p n + m. The two solutions of these 
equations of map [5] are 



1. 



-1 

2mn 



m 1 ■ <P°~ 
, p - — arcsm — 
2n \ y 



(27) 



and/3 (2) = l/2-/3 (1) , for each integer m = ±1, ±2, ±3, . . ., 
and 



2. 



2mn 



ff» = \ - ^arcsin(^) , (28) 



and/? (2) = 3/2 - = 1 + (1/2tt) arcsin(0cr/y), for each 
integer m = ±1, ±2, ±3, . . .. and arcsin is evaluated on 
the fundamental domain [-n/2, n/2] where it is uniquely 
defined. 

For positive and m = ±1 in Eq. [27] we find that corre- 
sponds to the hyperbolic (saddle) fixed point, whilst 0^ cor- 
responds to the attracting fixed point, as is illustrated in Fig. 
0(a). For the negative cp the same conclusion applies, respec- 
tively. For the whole sequence of larger | m |> 2 we have a 
sequence of fixed points but they can not be seen. The sad- 
dle fixed points given by Eq. (I27H28I) are marked by a star and 
are shown in the upper part of Fig. |4]a) while the Eq. (|28T> 
gives the saddle fixed point shown in Fig. |4ja) in the lower part, 
also marked by a star. The two branches of the unstable man- 
ifolds of each saddle are are obtained via iteration of the map 
P with appropriate initial conditions. The upper [Eq. (|271 »1 and 
lower [Eq. (l28t l unstable branches are marked by u and shown 
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Figure 5: (Color online) (a) Attracting fixed points for m = -1 represented by a upperward red star, downward green star and a chaotic attractor shown in black, (b) 
Their corresponding basin of attraction. The control parameters used in (a) and (b) were y = 10~ 2 and <x = 0.15. The region in black is the basin of attraction of the 
chaotic attractor and the region in red color is the basin for the upper attracting fixed point and the region in green corresponds to the basin of attraction of the lower 
attracting fixed point. 



in red in Fig. [4] They converge to the attracting fixed points, 
while the downward [Eq.(l27l)1 and upperward [Eq. (|28T »1 branch 
generate the chaotic attractor and are shown in green. On the 
other hand, the stable manifolds consist basically in trajecto- 
ries heading towards the saddle points and the construction of 
the stable manifolds is slightly different from those of the un- 
stable manifolds where we have to known the inverse of the 
mapping P, which we denote by P~ l . The basic procedure is 
that P~ l ((j) n+ \,p n+ i) = (<p n ,fin), thus the map P~ l is given by 

I. _ <pn + i-y sin(2flft„) 
B-B + - ' • (29) 

Once we have found P~ l , the stable manifolds, which corre- 
sponds to the borders between the basins of attraction for the at- 
tracting fixed points and the chaotic attractor are shown in blue 
and violet in Fig. [4] By reducing the value of the dissipation 
parameter cr the unstable manifold of each saddle point touches 
simultaneously the stable manifold and as a consequence the 
chaotic attractor as well as its basin of attraction are immedi- 
ately destroyed. Such an event is called as a boundary crisis 
JIM EES Eh and happens simultaneously for the manifolds 
arising from both Eqs. (l27l> and (1281 ) as shown in Fig. HJb) 
and is characterized for the first time in this model. The control 
parameters used in Fig. |4]were y = 10" 2 and: (a) cr = 0.15 (im- 
mediately before the crisis); (b) cr = 0.143 (immediately after 
the crisis). 

Before the crisis and for the combination of control param- 
eters used in Fig. |2a) (y = 10" 2 and cr = 0.15), there are 
three different attractors, namely: (i) two attracting fixed points 
related to m = ±1 (the central regions where the red spirals 



are converging to) and (ii) a chaotic attractor [see Fig. 0(a)]. 
One might expect that there must be three different basins of 
attraction. This is indeed true as one can see in Fig. 0b) . The 
procedure used to obtain the basins of attraction for the chaotic 
and attracting fixed points consist basically in iterating a grid of 
initial conditions in the plane j3 x and look at their asymptotic 
behavior. We have used a range for the initial ft as ft e [0, 1] 
and e [-0.15, 0.15]. Both ranges of j3 and were divided in 
1000 parts each, leading to a total of 10 6 different initial condi- 
tions and for the values of the control parameters that we have 
considered, each combination of (fi, 0) was iterated up to 10 5 
times. After the crisis, the time evolution of an initial condition 
in the corresponding region of the chaotic attractor before the 
crisis leads the orbit to make an incursion towards the region 
of the, now, chaotic transient until being captured by one of the 
two sinks. The structure of the two basins of attraction belong- 
ing to the two attracting fixed points is very complex and it is 
shown in Fig. so that it is difficult to decide which initial con- 
dition will go to which of the two attractors, as there are some 
obviously some regions in the phase spcae where both basins of 
attraction are dense. 

We may therefore argue that an initial condition chosen in 
the region of the phase space that produced a chaotic dynam- 
ics before the crisis, now gives rise to a chaotic transient. We 
denote as transient the corresponding number of iterations that 
the orbit typically spends until it finds the appropriate route to 
one of the attracting fixed points. The transient is described by 
a power law of the type 

n t oc^, (30) 
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tive Fermi-Ulam model [49] . The Fermi-Ulam model consists 
of a classical particle that is confined to bounce between two 
rigid walls. One of them is fixed and the other one is period- 
ically time- varying. Despite the remarkable difference of the 
two models (the main difference being the power of the term in 
the denominator of the two mappings, which is 1 in the Fermi- 
Ulam model and 1/2 in our model), the chaotic transient af- 
ter the boundary crisis marking the converging to the attracting 
fixed points follows the same rule in both models and is given 
by the same exponent = -2. 
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Figure 6: (Color online) Basins of attraction of the two sinks observed in the 
phase after the boundary crisis. The interwoven structure is evident. The control 
parameters used were y = 10~ 2 and <x = 0.143. 
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Figure 7: (Color online) Plot of n t vs.// for a fixed y = 0.01. A power law fitting 
gives £ = -2.09(7). 



3. Conclusions 

We have studied the problem of a charged particle in the 
electric field of a wave packet. We have confirmed through an- 
alytical arguments that close to the invariant spanning curves 
the dynamics of our model can be described by the standard 
map. Such approach allows us to find the position of the first in- 
variant spanning curve (chaos border) as a function of the con- 
trol parameter y. Once we have found that the position of the 
first invariant torus changes as a power of the control param- 
eter y with the exponent 2/3, we have studied the behaviour 
of the chaotic sea close to the transition from integrability to 
nonintegrability (small values of the control parameter y) using 
scaling arguments. We have shown that there exists an analyt- 
ical relationship between the critical exponents, namely, accel- 
eration exponent, saturation exponent and crossover exponent. 
Our scaling hypotheses have been confirmed by the perfect col- 
lapse of all curves onto a single universal plot. Finally, for the 
dissipative dynamics we have found the expressions for the sad- 
dle fixed points and constructed the corresponding unstable and 
stable manifolds. We have shown that increasing the dissipa- 
tion, the unstable manifold touches the stable manifold and the 
chaotic attractor as well as its basin of attraction are completely 
destroyed. We have shown that such a destruction is caused by 
a boundary crisis. Finally the transient marking the approach to 
the sink is characterized by a power law with exponent = -2. 



where \i - cr-cr c with cr c < cr. For a fixed y = 0.01, our numer- 
ical simulations lead to the critical value of cr c = 0.143056 . . .. 
Figure |7] shows the behaviour of the average transient plotted 
as a function of \i considering a fixed y - 0.01. A power law 
fitting gives the exponent £ = -2.09(7). The procedure used 
to simulate the transient consists in evolving in time an ensem- 
ble of B = 5 x 10 3 different initial conditions randomly chosen 
along the region of the chaotic attractor as observed just before 
the crisis, and calculated the time n\, where i is a member of the 
ensemble, that the particle spends until reaches one of the two 
attracting fixed points. The average transient is calculated as 

1 B 

i=\ 

where i specifies an orbit of the ensemble. The exponent £ = 
-2.09(7) has, within an uncertainty error, the same value as the 
exponent obtained for the boundary crisis observed in a dissipa- 
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